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1. Introduction 

Quantum entanglement [U El El H] is an important aspect of quantum mechanics that tells 
us about the quantum correlation of two particles or subsystems spatially apart. When 
a composite quantum system is in a pure non-product state, then even if the subsystems 
are spatially far apart and non-interacting, the measurement on one subsystem affects 
that on the other instantaneously. This "spooky action at a distance" later gave birth to 
the term "entanglement" . This phenomenon was first marked by Einstein, Podolsky and 
Rosen in a gedanken experiment [5] , known as the EPR paradox. In their paper, they 
considered two particles which interacted for some time and showed that it is possible 
to measure the conjugate non-commutating quantities, like position and momentum, 
simultaneously, in violation of quantum mechanics. Later on, it was resolved and 
resulted the idea of the quantum entanglement which indicates the presence of inherent 
quantum effects, not just correlations between the two particles. In fact, quantum 
entanglement is a resource for quantum computation with no classical counter-part. 
Here, in this paper, we consider two interacting particles in a pure state, like an EPR 
pair, or more specifically in the two-particle ground state, in which case we can expect 
to observe the entanglement. 

A quantum phase transition (QPT) occurs at zero temperature and at the QPT the 
ground state energy is non-analytic with respect to some parameter in the Hamiltonian. 
A QPT is fully governed by quantum fluctuations and hence one would expect that the 
quantum entanglement would show special signatures at the QPT. Even it is found that 
the entanglement entropy behaves in different ways for a first order and a continuous 
transition. The critical behaviour of the entanglement entropy is drawing much attention 
now-a-days[6] and has been investigated for different spin models [H El [7j as well as in 
continuum systems [8]. 

To quantify the entanglement, various definitions of entanglement entropy are used, 
among which the von Neumann entropy [1] El [7] is the most common. Recently it was 
shown that the von Neumann entropy of two particles has a dlnn behaviour at the 
quantum critical point (QCP) of unbinding in dimensions 1 < d < 4. This has been 
established analytically for a 3D potential well[9j. Here the QCP is attained when the 
inverse length scale k, the inverse of the width of the wave function, approaches zero. 
This is achieved by tuning the potential or the mass. Also it is shown that this divergence 
is essential for the criticality and linked to the reunion behaviour of two polymers in 
the equivalent classical statistical mechanical problem of polymers. In this paper, we 
study how the von Neumann entropy for a long-range potential, like 1/r 2 , changes as 
one varies its strength and sign. The equivalent classical statistical mechanical problem 
involves two directed polymers interacting at the same contour length like a DNA with 
native base pairing but with an additional 1/r 2 interaction. The polymer model has 
been studied using renormalization group in Ref . pTOj [TTj . 

This paper considers two particles interacting through a three-dimensional inverse 
square law potential, and finds the quantum entanglement between the particles. Here 
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we use particle partitioning]!)]. The Hamiltonian for the two particles we shall be using 
is, 

where mi, and are the mass, the position and the momentum of the zth particle, 

( -V Q , for r < a, (V > 0), 
V{r) = { 2^g t (2) 

{ fOT r > °> 

is a central potential, and fi = m^/ (mi +777,2) is the reduced mass of the two particles. 
We take = 1. 

The inverse square potential is of immense importance in quantum mechanics [12]. 
It is at the boundary of the short and the long range potentials. For potentials decaying 
like r~ p in three dimensions, there is no finite bound state if p > 2, while for a slower 
divergence, i.e., p < 2, there is a finite negative lower bound in energy. For an attractive 
potential —g/r 2 (g > 0), the kinetic and the potential energies are of the same order near 
small r and so the bound state spectrum depends on the value of g. A manifestation 
of the borderline case is in the scale-free nature, H(Xr) = \~ 2 H(r). This makes, g, 
the dimensionless strength of the potential, a "marginal" parameter in the RG sense in 
all dimensions. The singularity of g/r 2 at the origin prevents discrete bound states. A 
suitable modification of the potential at small r, e.g. by putting a cut off and replacing 
the potential by a short range finite attraction near origin, gives discrete bound states. 
This is done in Eq. [2j 

It is established in quantum mechanics that there is no finite energy ground state 
for g > 1/4. For g < 1/4, the wave function is normalizable and the bound state 
energy can be obtained by the standard procedure. In this regime of g, the unbinding 
transition is induced by tuning the strength of the short-range potential near r = 0, 
depicting the quantum phase transition. The unbinding transition in this long-range 
interaction is a unique example of a QPT whose type can be first order (g < —3/4), 
critical but non-universal (g > —3/4), and even Kosterlitz-Thouless type (g = 1/4) [11]. 
The solvability and the wide repertoire of QPT behaviour make this model an ideal 
terrain for exploration of the nature of entanglement entropy around a This is 

what we set to do in this paper. 

A phase transition is defined as a singularity in the energy, associated with diverging 
length scales. In this sense the quantum unbinding transition is a genuine phase 
transition. This QPT exists because time of infinite extent plays a role in quantum 
mechanics. It becomes clear in the path integral formulation. The quantum problem 

I What makes the unbinding transition special is the change in the topology of the underlying space. 
Note that the wave function \ifj) belongs to a space Hb — > H c ® Hi where H c is the Hilbert space for a 
free particle while Hi is the separable space of square integrable functions. This space for the bound 
state goes over to the space Hjj = H C ®H C for two particles (scattering states) . Another way of seeing 
this is via currents. The space Hb allows only one current for the CM while Hjj allows two currents 
for the two unbound particles. In this paper we stay in the Hb space on the bound side. 
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-2 


-3/4 





1/4 


> 1/4 


A 


1.5 


1 


1/2 





imaginary 



Table 1. Conversion table for g-X. Real values of A occur for g << 1/4. 



can be mapped onto an equivalent classical statistical mechanical problem of polymers 
under the imaginary time transformation (it — > N). The time in the quantum problem 
then becomes the length of the polymer, N, Green's function maps on to the partition 
function and the ground state energy corresponds to the free energy per unit length. The 
interaction between the polymers means the interaction of a pair of monomers at the 
same index along the length of the polymers as in DNA base pairing. This is equivalent 
to the same time interaction of two quantum particles. The equivalent classical problem 
in the context of melting transition of two polymers interacting via a potential like Eq. [2] 
has been discussed in Ref. pi)] which reveals that the results of the quantum problems 
can be recovered from such studies. Like the quantum particle making excursion inside 
and outside of the well, the polymers also come closer, they reunite, and move further, 
forming swollen bubbles. The entropy of a bubble of length iV is 

In ft (iV) = Na -tflniV, (3) 

where Q(N) is the reunion partition function of two polymers starting together, reuniting 
anywhere in space again at length N, <To is the bubble entropy per unit length and \I/ is 
the reunion exponent. The details can be found in Refs. [9j QUI [13] . 

The binding-unbinding transition of polymers, for DNA melting, has been studied in 
the context of the necklace model of polymers and it is found that the reunion exponent 
\P determines the order of transition [13J. The phase transition occurs if ^ > 1. The 
transition is continuous if 1 < ^ < 2, while it is first order for ^ > 2. In three 
dimensions, the reunion exponent is given by|10] 



#=1 + A, With \ = J^-g, (4) 

where the dependence on g, a bit counter-intuitive, is a consequence of its marginality. 
Here also we use the parameter A(> 0) because of its occurrence in the sequel. Table 1 
gives the correspondence between g and A for easy reference. 

The phase diagram and the lines of RG fixed points are shown in Fig. [TJ This plot 
shows the phases in the g-u plane, where 

u = —Vq a 2 , (5) 

in the unit of 2fi/h 2 = 1, is the dimensionless short range potential in which the two 
particle state is in. The fixed points shown here are obtained from the renormalization 
group analysis done in Ref. [10]. The red line for u < —0.5 shows the unstable fixed 
points across which the unbinding transition takes place, and the black continuous line 
for u > —0.5 shows the phases by stable fixed points. For g < —3/4, the bound-unbound 
transition is first order for ^ > 2, which is indicated by the red continuous line ending 
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Figure 1. g vs. u phase diagram. The plot shows the phases and the RG fixed points 
in the g-u plane (u = — Vq a 2 ). The red curve below g = 1/4 and u = —0.5 show the 
binding-unbinding transitions governed by a line of unstable real fixed points. The 
transition is first order for g < —3/4 and second order for —3/4 < g < 1/4. This 
line is the transition line in the limit of zero range potential (a — > 0, Vq — > oo, with 
u =constant) The black continuous curve for u > —0.5 shows the locus of stable fixed 
points representing the unbound phase. The dashed-dotted line at g = 1/4 is the 
boundary beyond which the fixed points are complex. 

at the symbol * at g — —3/4, or, A = 1, a multicritical point. After that, the transition 
is continuous upto g = 1/4 with \& < 2. Beyond that, where A is imaginary, there is 
no real fixed point, and the system is in a bound state. Across the g = 1/4 line, with 
u > —0.5, a Kosterlitz-Thouless type phase transition from the bound to the unbound 
state can be induced by tuning g. The two regimes, ^ < 2 and ^ > 2 are governed by 
different behaviours, with additional log-corrections at \1/ = 2. 

We find that the entanglement entropy also carries this signature of the speciality 
of g = —3/4 or A = 1. The entropy in the three different regimes, A < 1, A = 1 and 
A > 1, scale in different manners. We establish that A = 1 behaves like a multicritical 
point for the entanglement entropy too, controlling both the first order and the critical 
behaviour in the whole range —2 < g < 1/4. 

The outline of the paper is as follows. In section 2, we describe our model and the 
method by which we calculate the von Neumann entropy. The analytical results are 
presented in section 3 and the von Neumann entropy is calculated for A < 1. Next we 
present the exact numerical results done in mathematica and discuss the behaviour 
of the entropy and its scaling in section 4. Finally we conclude in section 5. 

2. Model and method 

Eq. |2] is used for our study. The detailed nature of the short-range potential is not 
important and we take it as a simple square well potential. We concentrate in the range 
< A < 3/2. 



Entanglement entropy for inverse square interaction 6 

The von Neumann entropy (S) is denned as 

S = -Tr(p\np) (6) 
where p is the reduced density matrix for the ground state \ip) of a two particle system, 

p(n,rl) =Tr 2f? (l,2) = / d 3 r 2 (n, v 2 \^)W 1 , r 2 >, (7) 



obtained from the two particle density matrix g(l, 2) =| if))fy> | by integrating out particle 
2. The density matrix in hip is to be made dimensionless by appropriate powers of cut 
off a. 

As the interaction is translationally invariant, the reduced density matrix p(r, r') = 
p(r — r') satisfies the eigenvalue equation[9], 

J d 3 r' p(r — r') exp(— zq ■ r') = p(q) exp(— zq • r). (8) 

The eigenfunction is exp(— iq • r) with the eigenvalue 

P(q) = l0(q)| 2 , (9) 

where we have put the center of mass wave vector to be zero without loss of generality. 
Here 0(q) is the normalized momentum-space wave function of the relative co-ordinate. 
Therefore, the von Neumann entropy reads, 

S = - J d 3 q|0(q)| 2 ln|0(q)| 2 . (10) 

The reduced density matrix in the basis of momentum states |k) has the form 

p = Jd 3 k |0(k)| 2 |k)(k| ^ J d 3 k |k)(k|, (11) 

which makes the mixed state characteristic explicit. Eq. [TT] allows us to define p as a 
thermal density matrix with an entanglement Hamiltonian H ent at a fictitious inverse 
temperature with Z as the partition function. This thermal correspondence makes 
the von Neumann entropy equivalent to the Gibbs entropy of H ent . 

In Eq. [HI H ent is a c-number. Consider the canonical partition function of a free 
particle at temperature T, Z ~ J d d qexp(—(3H) ~T d / 2 , where H = h 2 q 2 /2m. Then 
the entropy becomes, S = \nZ ~ InT, which for very low temperature, T — > becomes 
negative. In another way, one gets a constant specific heat C from the equipartition 
theorem, which then gives a logarithmic dependence on temperature of the entropy, 
S = J (C/T)dT ~ InT. The Sackur- Tetrode constant, the dimensionless entropy of 
one mole of an ideal gas at temperature T = IK and one atmospheric pressure, is a 
fundamental constant[15]. Its value is —1.1648708. Note that this fundamental entropy 
is negative. Classical harmonic oscillator is no exception. It is well-known that the 
condition S > does not hold for the classical continuous statistical mechanics[16j. 
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3. Analytical results 

In this section, we derive the asymptotic behaviour of 4>(q). In particular we find that 
the entropy is dominated by the outer part, i.e. the excursion in the classically forbidden 
region, if the unbinding transition is critical. This happens for < A < 1. For first 
order transition, the inner part also contributes significantly. 

To find out the von Neumann entropy of two interacting particles in the ground 
state, we first write the wave function as a product of that for the center of mass (CM) 
and for the relative coordinate: 

^(r 1 ,r 2 ) = $(R)^(r), (12) 

where $ and tp are the wave function in the CM and relative coordinates respectively, 
defined by, 

_ miti + 1712^2 

R= and r = ri — r 2 . (13) 

nil + n^2 

We proceed with the radial part putting the CM momentum as zero. The CM is 
completely delocalized in space. 

The ground state has zero angular momentum and therefore, only the radial part 
of <p(r), R(r), in the spherical polar coordinate is needed. For this s-state, the radial 
part of the Schrodinger equation then reads [13]: 

Q2 p n an 

7nr + -7T + (V + E)R = 0, for r<a, (14) 

d 2 R 2 dR /a \ 
and, ^— + -— + (4 + ^)^ = 0, for r>a, (15) 



where E is the ground state energy of the particle describing the behaviour of the two 
particles in relative coordinate. The radial part of the wave functions in the relative 
coordinate is then obtained by solving Eqs. H3J and \T5\ 



— sinfcr, for r<a (16a) 

r 

—j=H^ \inr), for r > a, (16b) 



R{r) 



where A and B are the normalization constants, H^' is the Hankel function of first 
kind, and 

k 2 = V - \E\, k 2 =\E\. (17) 

We choose A to be positive and it is given by Eq. HJ In the limit of k — > 0, the unbinding 
transition takes place. This makes our interest in studying the von Neumann entropy 
in this limit. 

At r = a, the continuity of the wave functions gives, 
A B 

— sin/ca = —=H^(iKa), (18) 
a \ a 
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and the matching condition of the derivative of the wave function gives, 

ifR(zKa) 1 . 

akcotak = ina — j-, AH — , (19) 

H^(iKa) 2 K J 

which determines the value of k for a given k. Given the values of A and a, one can get 
the threshold or minimum value of k, k m , for just one bound state. For k = 0, 

ak c cot ak c = — — A, (20) 



is the condition for the transition point when the ground state energy E — > 0. Eq. 
has always a solution for A > 0. 

Now consider a small deviation from the critical value of k, k = k c — 5 where 
6 ~ Vo — V c . Then, from Eq, 

f (na) 2X , for A < 1, 
(ak c — aS) cot(ak c — a5) ~ < (21) 

I (na) , for A > 1, 

or, 



„ { 5 1/x , for < A < 1, 

\E\ = K*~\ (22) 
' ' {5 + 0(5 l ^ x -V), for A>1. V ; 

These show that as V — > V c = k%, E remains continuous, as it should. For A < 1, E 
approaches zero tangentially, while for A > 1, there is a nonzero slope at k — 0. This 
discontinuity of slope classifies the A > 1 transition as first order. Despite that, the 
higher derivatives on the bound side d n E/d5 n would show divergences like a critical 
point. 

The normalization constants A and B are found by using the continuity condition 
and taking the limit k — > (see Appendix for details), 

_ f (aK) 2 - 2X /a, for A < 1, 

I 1/a, for A > 1, 

and 

„ [ k 2 , for A < 1, 

LB 2 ~< (24) 

1 ' \ K 2 {aK) 2X -\ for A > 1. 1 J 

At A = 1, there are log corrections which we do not get into here. The log correction 
appears in the Necklace model for polymers whenever the reunion exponent ^ (Eq. |4]) 
is an integer. The log appears in Eq. [19] via Hq for A = 1. Now one knows the full 
wave function and its limiting k behaviour. 

The reduced density matrix has eigenvalues |0(q) | 2 , where q is the momentum space 
variable. To get these eigenvalues, the Fourier transformation of the wave-function needs 
to be done, 

m = -j^f 2 f d 3 r e^R(r) = <k(q) + o (g), (25) 
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where the subscript i,o refer to the inner (r < a) and the outer (r > a) part. The 
Fourier transform of the inner part (Eq. I16aj) is 



Hi) 



A 1 
q V2tt 



sin(fc — q)a sin(g + k)a 



k 



q + k 



(26) 



and of the outer part (Eq. I16bl) is 

4>o{q) 



m «-</> £ r + * 

4 2 



7T 



2^1 



A 

2' 



A 3 
2' 2' 



K 2 



dry?" /vx(/cr) 



(27) 



'o 9 

where 2 -Pi is the hypergeometric function and K\ is the modified Bessel function. The 
last integral in Eq. [27] is convergent for all A < 1 and therefore can be ignored in the 
k — > limit. 

The limiting small k dependence of the inner and the outer parts of the wave 
function from Eqs. [26] and [27] are of the form 

,l-A 7 3/2. 

(28) 



Uq) 




l fi(aq), 



fi(aq), 



if A < 1, 
if A > 1, 



and 



K 



K 



(29) 



/a(?/k), if A < 1, 

(a K ) A -V A (g/ K ), if A>1, 

where fi and fx are well-behaved functions. Eq. [28] is for large q. 

From Eqs. [28] and [29], we see that the double limit k — > 0, A — > 1 is singular because 
of the term (an) 1 ~ x . The dependence of entropy on the order of the limits is discussed 
later on. This identifies (n = 0, A = 1) as a special point. From this we also identify 
(1 — A) ln(an) as an appropriate scaling variable. This scaling variable will occur below 
in the analysis of the numerical results. 



For A < 1, i.e., 1 - A > 0, 



1 an 



,1-A 



— > as k — > and therefore, the contribution of 



outer part dominates over the inner part in the von Neumann entropy. Without much 
loss, one can then write the entropy with the outer part only (Eq. 129]) . 



S 



d 3 q |0 
3 ha. an + C\. 



2 hi|0 o 
(for 



a< r 



with 



c x = / dxx 2 f x (x) lnfx(x). 



(30) 



(31) 



We introduced a in Eq. [30] to make the argument of In dimensionless, mentioned earlier. 
As per our interest, we extract the K-dependent term and call the rest C\, which is a 
function of other parameters. The main result is that there is a log divergence of S as 
k — y 0. 
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4. Exact numerical results 



To study the nature of the entanglement entropy, over the whole range of A we take 
recourse to exact numerical calculation using mathematica for the 3-dimensional 
potential well. We cross-check our prediction of Eq. [30] and then show a multicritical 
scaling that covers the range < A < 3/2. 



4-1. Protocol 

Although Vo is the tuning parameter, it is more convenient to use the length scale as the 
independent parameter. With this treading of k for Vo, our protocol is like this: Given 
the values of k and A, the value of k m was determined from Eqfl9| with k m < tt that 
assures us the ground state. As k — > 0, k m — > k c . Then the corresponding normalization 
coefficients A and B were found by using the normalization condition and the continuity 
equation, i.e., by doing the r-integrations of the inner and the outer parts of the wave 
function in Eq. IA.2I These constants are used in the Fourier transformed inner and outer 
parts of the wave function, Eqs. |26]and[271 to calculate the von Neumann entropy. In the 
final integration for S = — Tr plnp, we put an upper cut off making sure that the final 
numbers are independent of this choice of cut off. Also the intervals of the integration 
range have been chosen carefully especially for q ~ k. This gives numerically exact 
numbers for the entropy for the given k and A. This procedure is repeated for various 
A and k. We set a = 1. 



20 



oo -40 



-60 



-80 



-100 



L--A""i 

»•♦-< 


^.A-A-A-^ 
^-O -0_ O-< 


. • A- A • A- ^ 


..A--A-A--. 


■if 




* * K=10 A 
(-*.•-'-• K=10" 2 
— K=10 J 
— • K=10 4 




- H- .-+-+-H 








+ " + K=10" 5 
-• K=10" 6 
""- K-10" 7 






..+■■*■•■*'' 


.» 4 .' 
' .* 'Ill 

,"';/!! 
y /'I 


i 
i 
i 




- ' K=10" 8 

-9 

- ' K=10 










i 


- K=10- 10 

- k=io-'; 

K-10"' 2 
" K=10- 13 












i 












i 




- K=10 14 



0.4 0.6 0.8 1.0 1.2 1.4 

X 



-u 



Figure 2. S vs. A for various k. The plot shows that the entropy diverges for A < 
as k — V 0. The dashed line marked as k — > is the expected behaviour of the entropy 
for A > 1. 
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4-2. Behaviour of the von Neumann entropy S 

4-2.1. A dependence: The plots of the numerical values of the von Neumann entropy 
S against hiK and A show different behaviours of entropy in different ranges of A, viz., 
A < 1, A > 1 and A = 1. 

Let us first look at the plot of S vs. A in fig. [2j where different lines represent 
different values of k. For A < 1, the von Neumann entropy for small k saturates to a 
negative value as A is varied and that saturation value depends on the value of k. The 
smaller the value of k, the more negative is the entropy, and k — > takes the saturation 
value to negative infinity. The long-range part of the potential is attractive for A > 0.5 
and repulsive otherwise. But the entropy shows no signature as it crosses A = 0.5. On 
the other hand, for A > 1 where the transition becomes first order, the entropy does not 
decrease much with k, rather becomes independent of k. It remains finite for A > 1 and 
diverges at A = 1 like the black dashed curve in Fig. [2j S becomes positive at A ~ 1.3. 
It seems that this point has no significance otherwise. 

4-2.2. k dependence: The behaviour of the von Neumann entropy with A and k becomes 
more clear when one looks at the plot of S vs. k (Fig. [3]). This plot shows the different 
characteristic behaviours of S in the three distinct ranges of A: A < 1, A = 1 and A > 1. 
For small k, all A < 1 curves have slope 3 when plotted against ln/t, i.e. for A < 1, the 
entropy is of the expected form 3 In k + c\ which is shown from analytical calculations. 
To get 3 In k, one has to see below some value of k, and as A approaches one, even 
smaller k needs to be considered. But no matter how close to 1 is the value of A, one 
gets 31n/t until A < 1. Exactly at A = 1, the slope changes suddenly to 3/2 and hence 




In k 



Figure 3. S vs. In k for various A. For comparison, 3 In k and | In n are shown by black 
lines with the symbols + and x . The inset shows that the entropy is K-independent 
for A > 1. 
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A somewhat different behaviour is seen for the rest with A > 1 (inset of Fig. [3]). 
For small k, the curves reach a A-dependent constant value and do not change with k. 
Clearly the entropy has no K-dependence for A > 1 and it is finite. By definition, these 
constant values are c\ and S(\ > 1) = c\. So, we see that there are three classes: 



S 



3 In k + C\ 
- In k + ci 



for A < 1, 
for A = 1, 
for A > 1. 



(33) 
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Figure 4. The plot of c\ vs. A, showing a divergence at A = 1. 

4-2.3. On C\: Now we have knowledge of the k dependent part in the von Neumann 
entropy for different A. The next question is how the c\ behaves with A, and if they have 
different nature in different regimes of A. So, we collect the c\s according to Eq. [331 
except for A = 1, and plot against A. This plot (Fig. HJ) shows a divergence at A = 1 
indicating that (1 — A) is an important quantity. The data points are fit into the function 

CA = m + n/(l - A), (34) 

via m and n, and the fitted set of parameters are (4.52, 1.38) and (3.76, 1.48) for A 
greater and less than one respectively. The divergence of c\ at A = 1 leads to the 
possibility of a reduction of the slope of S from 3 to 3/2 when plotted against In k. 

4-2.4- Data collapse: We noted that for A > 1, C\ and hence the entropy itself, has a 
(1 — A) dependence and for A < 1, the entropy has a ln/t term with c\ = /(l — A). It 
was pointed out in Sec 3, below Eq. [29] that (1 — A) Inn seems to be a scaling variable. 
We therefore look at the plot of the entropy vs. (1 — A) ln/t. The entropy has different 
behaviours on the two sides of the A = 1 making it a special point. Also, it has a 
separate scaling behaviour. This drives us to plot (S\ — Si) / (| In k) against (1 — A) In k. 
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Figure 5. Data collapse: (5a — 5i)/|ln«; vs (1 — A)1iik. 



We see a good data collapse (see Fig. EJ) for various sets of data of Fig [2j Hence, one 
can write the scaling form of von Neumann entropy: 



(S. 



x 



St)/ -Inn 



A) In k). 



(35) 
■1 for 



Fig. [5] shows that (S\ — 5i)/|lnK reaches +1 for small enough k for A < 1 and 
A > 1. Once we get the scaling behaviour of the entropy at A = 1, the same away from 
this special point can also be obtained. 



5. Discussion and conclusion 



In this paper we studied the von Neumann entropy S, the most common measure of the 
entanglement entropy, for an inverse square potential in three dimensions. The entropy 
behaves in different ways for three different ranges of modified interaction strength A 
and given by Eq. [33 For A < 1 and A = 1, the entropy has a diverging nature as one 
approaches quantum critical point by tuning k. The behaviour of entropy is completely 
different for A > 1, where the ^-dependence of the entropy vanishes. There is a 
divergence in the entropy. These three distinct classes collapse onto a single curve when 
(S\ — Si)/§ In k is plotted against (1 — A) In k. This data collapse indicates that there is 
a common scaling behaviour of the entropy for any A and that A = 1 is special. Because 
of the diverging factor dependent on (1 — A), one has to be careful in taking required 
limit of k — y as that would give a log correction in entropy for A = 1. For g > 1/4, A is 
imaginary which we do not consider here. In this paper we focused on the multicritical 
point at A = 1. There is one more multicritical point at A = with KT transition. 
This case will be discussed elsewhere. It would be interesting to study the entanglement 
behaviour for a discretized version of the model. 

The nature of the divergence of the entanglement entropy at A = 1 depends on the 
path of approaching A = 1 in a A-k plane. Diagrammatically it has been shown in Fig. |6j 
If we take the limit k — y first and then A = 1, the entropy diverges like 1/(1 — A) 
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Figure 6. The path dependence of entropy. Two different limits of approaching 
{A = 1, K —> 0}. (a) First k — > and then A = 1 (red line). The entropy diverges like 
1/(1 — A), (b) First A — > 1 and then n —> (blue line). The entropy diverges like InK 
Thick line along x-axis for A < 1 denotes divergent entropy, (c) For A < 1, taking the 
limit k — y (black vertical line with arrow) leads to divergent entropy and S remains 
so along the horizontal stretch. 



(shown by red line (a) in Fig. [6]), and like Inn for the other way around (see the blue 
line (b) in the same figure). For A < 1, the k — > line corresponds to S = — oo, but 
for A > 1 the same line gives a finite value for entropy. The path dependence of Fig. [6] 
summarizes the features of the entanglement entropy with A = 1, /t = as a special 
point controlling the behaviour in its neighborhood. The data collapse of Fig. [5j then, 
suggests that the paths should be classified by the constant value of X = (1 — A) In k. 

For A < 1, restricting to the critical case, we see p(q) ~ |0(g)| 2 . These are the 
eigenvalues of the density matrix. Now the reduced density matrix p describes a mixed 
state, though the full ground state is pure. Being a mixed state, we may represent p as a 
"thermal" density matrix, p ~ exp(—j3H ent ), as done in Eq. [TTJ Since the entanglement 
spectrum is known, we have 

(5H ent w In | 2j F\| 2 w ~4, for q -> 0, (36) 

identifying (3 = 1/k 2 and H ent = q 2 /2. As mentioned before, for this H ent the Gibbs 
entropy is ~ | InT. Since in this case T ~ k 2 , we find the von Neumann entropy 
S ~ dlnK. 



Appendix A. Calculation of the normalization constants A and B 



The normalization constants A and B are found by using the continuity condition and 
taking limit k — > in the normalization condition, 



47T 



I A I sin kr dr + 



r\BH^\inr) | 2 dr 



(A.l) 
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Replacing A by B, by using the continuity equation, Eq. (fig]) , we get, 



7T 



27ra — — sin2a£; 
k 



a\H^\iKa)\ 2 



+ 47T 



r|ff A (zkx)| 2 dr 



sin 2 ka 

Now we use the form of the Hankel function in the limit k — > 

2 A r 2 (A) 



fff } (kr) | 2 ~ 



/tr 



-2A 



7r z 



|S| a = l. (A.2) 



(A.3) 



and rewrite the outer part integral in the normalization condition in a simpler form, 



mz. , 

POD 

/ r | H^\ir) | 2 dr 

J na 



Ka . 

2 



\Ht\ir, 



|2 2^r 2 (A) r _ 2A 



7T 



dr + 



x 2 A r 2 (A) 2Aj 

r~ 2A dr + 

7T 2 



r|iJ^(ir)| 2 dr 



r 2 (A) i 



I an 



,2(1-A) 



+ 



7T 2 2 X - A 1 - A 
Putting Eqs. (IA.4|) and flA.3|> in Eq. HA.2|) one gets 

'A 2 + a 2 fc 2 -l/4 



LB 



7T 



2 A +!r 2 (A) 



fc 2 



-2A 



1 1 - (a/c)^ 1 -*) 

^2 (T^A) 



(A.4) 



(A.5) 



which in the extreme limit of k — > gives the k dependence of B, 

„ [ k 2 , for A < 1, 

l^| 2 „ ' . 9 (A.6) 
' ' \ K 2 (aK) 2X - 2 , for A > 1. V ; 

Once B is obtained, the ^-dependence of the other constants A can be found by using 
Eq. (TT8]) as, 

. (a/t) 2 ~ 2A /a, for A < 1, 
\A\ 2 ~\ / 7 ' ' (A.7) 

1/a, for A > 1. 
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